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CN ■ Abstract 

^ | The symmetry approach is used for classification of integrable isotropic vector 

Q ' Volterra lattices on the sphere. The list of integrable lattices consists mainly of new 

equations. Their symplectic structure and associated PDE of vector NLS-type are 
discussed. 

^ ; 1 Introduction 

^ ■ We call vector Volterra lattices the equations of the general form 

V niX = f n V n+ i + g n V n + h n V n -i, neZ 

— ^ ■ 

> 

where V n are vectors and f n ,g n , h n are scalar functions depending on V n+ x, V n , V n ~\- The 
integrability is understood as existence of higher symmetries, that is the equations which 
are consistent with ([1]), but involve the larger number of neighbor vectors (preserving the 
same quasi-linear structure). The precise definitions are given in next Section. The goal of 
this paper is to classify integrable cases under the following assumptions: 

(i) the lattice and its symmetries are isotropic and shift invariant, that is their coeffi- 
cients depend only on the scalar products v m , n := (V m , V n ) = (V n , V m ) and this dependence 
is the same at each node; 

(ii) the lattice must be integrable independently on the dimension of the vector space 
and the nature of scalar product; 

(iii) all V n are of unit length, v n , n = 1. 

The shift invariance allows the use of shorthand notation with the discrete variable n 
omitted from subscripts, so that equation (CQ) takes the form 

V x = fV 1 + gV + hV- 1 (2) 

(subscripts x, t will be reserved for denoting derivatives, not shifts). Due to the other 
assumptions, functions /, g, h are related by equation 

vi,of + g + v Q -ih = (3) 



On 
(N 

(N 

o 



*L.D. Landau Institute for Theoretical Physics, la Semenov pr., 142432 Chernogolovka, Russia. 
E-mail: adlerQitp .ac.ru 



1 



and depend only on the scalar products t>i i0 , Uo,-i> which can be considered as indepen- 
dent variables. Therefore, the classification problem is reduced to finding of two functions 
of three variables, so that its complexity is comparable with the case of scalar Volterra 
lattices 

v x = f(v 1} v } v^i) (4) 

classified by Yamilov [I], see also the recent review article [2]. The whole method of solu- 
tion is also very close, since the necessary integrability conditions in both cases formally 
coincide (the difference is in the set of dynamical variables: t> m .n instead of v n ). In the 
continuous case, the general approach based on this remarkable observation has been de- 
veloped by Sokolov and Meshkov in the pioneering papers 0, H] devoted to the classification 
of KdV-type vector equations (including the anisotropic ones) on the sphere. Important 
classification results for some other types of vector PDE were obtained in [SlElEj, however 
the approach in these papers relied essentially on the polynomial or rational structure of 
equations. 

In principle, the classification problem for the lattices (j2J) can be solved without the 
unitary condition (iii). This constraint does not define an independent class of equations, 
but only a special reduction of the general problem. Indeed, it can be resolved by use of 
the stereographic projection 

_ 1-(U,U) 2 tt 

v -i + {u,u) eo+ i + (u,u) u > 

where eo is some fixed unit vector and U belongs to its orthogonal subspace. Vector U 
satisfies, in virtue of equation ([2]), some isotropic lattice U x = fU\ + gU + hU-i. Since the 
dimension of the vector space is inessential in our considerations, we see that any lattice 
on the sphere corresponds under this mapping to some lattice in the free space. On the 
other hand, this lattice for U is not arbitrary: it must admit the reduction (U,U) = 1 since 
U = V under this constraint. This reduction brings back to the original lattice. 

The paper is organized as follows. Section [5] contains a concise explanation of symmetry 
approach and derivation of the sequence of integrability conditions in the form of conser- 
vation laws. These are used in Section [H] which is the main technical body of the text. 
All lattices (J2j) are divided there into two subclasses; the first one is analyzed thoroughly, 
while the second one is poor in answers and its presentation is more brief. The results of 
classification are presented in Section [U The rest of the paper contains some discussion of 
associated PDEs and symplectic structures. 



2 The necessary integrability conditions 

The symmetry approach to classification of integrable equations had been developed in 
80's, see e.g. [H |TU] as general sources, [IH] for a modern account on the discrete case 
and review articles [HI [2] for detailed references. The lattice §2§ is called integrable if it 
possesses an infinite hierarchy of the symmetries of the form 

V tk = P {Kk) V k + V^'^Vu-! + ■■■+ p^' l ~ k) V^ k + p( fe >- fe V_ fe (5) 

with coefficients depending on the scalar products of Vt, ... , VLfc. It is easy to see that 
the compatibility condition splits over the vector variables V n resulting in the commutator 
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relation 

D x (P^)-D tk (F) = [F,P] (6) 

for scalar operators 

F = fT + g + hT~\ P {k) = p ^T k + ■■■+ p^- k ^T- k 

where T denotes the shift operator n \— > n+1. This allows to use the necessary integrability 
conditions established in the scalar case (j4j) by Yamilov [1] , with operator F instead of the 
linearization operator /* = f Vl T + f v + f v _ 1 T~ l . For sake of completeness we repeat very 
briefly the derivation of these conditions. 

Equation is equivalent to a set of equations for the coefficients of P^ k \ One pair of 
equations defines explicitly the leading coefficients 

p {k ^ = / fc _x . . . frf, p^ = ah. k+1 . . . h^h, 

while solvability of the rest equations provides some sequence of necessary conditions to 
integrability of the lattice. These conditions do not depend actually on the order k of the 
symmetry. More precisely, let equation (jHJ) can be solved, at some k = K, with respect to 
21 coefficients p( fc . ±fc ) ; p( fc - ± ( fc - 1 )) ) . . . ; p(*>±( fc -0) ; where k — I > 1. Then it can be solved with 
respect to these 21 coefficients at any k > K. Moreover, the coefficients of one symmetry 
are expressed through the coefficients of the other one by explicit formulae. In order to 
prove this, it is sufficient to notice that the term D tk (F) in the l.h.s. of (EJ) affects the 
computation of the coefficients p^' 1 ^^ p( k >°\ p^^ 1 ) only, and that the special form of the 
leading coefficients written above allows to approximate by the formal power series 
(p( K )^j k / K _ This brings to the following statement. 

Statement 1. If the lattice (TJ|) possesses an infinite hierarchy of higher symmetries then 
the equations 

L x = [F,L), L — a^T + a® + a^T' 1 + a^T~ 2 . . . 
L x = [F,L), L — a^T- 1 + a® + a^T + a^T 2 . . . 

are solvable with respect to the coefficients aP> } aP) depending on v mM . 

The series L, L are called formal symmetries. In turn, the equations for their coefficients 
can be rewritten further as the sequence of conservation laws 

D,(pW) = (T-l)(aW), D x (p®) = (T~ l - j = 0,1,2,... (7) 

More precisely, if the lattice ([2D possesses the symmetry of order k, then equations (|7j) can 
be solved with respect to for j — 0, . . . , k — 2. The densities are expressed 

explicitly by certain recursive algorithm in terms of the lattice coefficients and previously 
found This algorithm relates p^> with the residue of V defined as the free term of 

power series in T (the formula res[v4, B] e Im(T — 1) can be proven). However, in practice 
we will need only few several conservation laws and the corresponding formulae can be 
derived straightforwardly. 
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Statement 2. Let the lattice (0|) be integrable, then equations are solvable for the 
following sequence of the densities p^\ p& : 

p®=\ogh, (8) 

P {1) =g + ^\ (9) 

P (2) = Ai + ^(p (1) ) 2 + ^ {1) . (10) 

Proof. The equations for the coefficients aC-^aWaM a< 2 > are: 

= M" 1) -/ia ( - 1) , 

aW = /ai 2) - f^ 2 a® + ^a^ - y_ ia « + ha® - ha®. 

The first equation implies a^~ x > = f, without loss of generality. Then the second equation 
takes the form (log/% = (T — l)(a® — g), so that we obtain the density p® and the 
formula for the next coefficient of the formal symmetry: a® = g + a®. Accordingly to the 
third equation, this coefficient may be taken as the density p^ and then = h+a^ / f-\. 
The last equation can be brought to the form 

(hU + \{P W ? + * {1) ) x = (T- l)(f^ 2 a® + atop®) 

after multiplication by /_i and taking into account the previous equations. The second set 
of the densities is obtained immediately due to the symmetry n —>■ —n. □ 

Remark 1. In addition to the higher symmetries, existence of the higher order conservation 
laws is another characteristic feature of integrable equations. It is possible to derive some 
integrability conditions from this property as well. This leads to the notion of formal 
conservation law 

S X + SF + F T S = 0, S = s® + s^T- 1 + s®T~ 2 + ... 

where (glP) t := T~^a and coefficients s^' depend on v m>n . Solvability of this equation is 
equivalent to the sequence of conditions of the form 

p^ = (T-l)(aW), j = 0,1,2,... (11) 

In particular, 

p® = log(-///i), p® = 2g + D x (a®). 

It can be proven that conservation laws (171) are equivalent in virtue of conditions ( TTTT) . that 
is p {j) + const p ij) e C © Im(T - 1). In some classification problems use of these additional 
integrability conditions may lead to a crucial simplification or even to a shorter list of 
equations. In particular, these conditions were used by Yamilov in his classification of the 
scalar lattices (jlj) (see footnote on p. 567 and Theorem 22 in [2]). It turns out, however, 
that in the vector case these conditions are of minimal value and it is possible to dispense 
with them (in all found lattices they are fulfilled automatically). 
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/'-= hf-^-ip^f + o- 
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Returning to the characteristic equation we notice that solvability of the first pair 
of integrability conditions ((7j), (jSJ) allows to find the coefficients p^ k > ±k \ p(. k >±(k-i)) f 
symmetry. At k = 2 this defines the symmetry completely, due to the constraint (V, V) = 1 
which implies 

V2,0P (2 ' 2) + V l>0 p^ + P i2fi) + V ,- lP ^- l) + ^ ,-2P (2 '- 2) = 0. 

The straightforward computation shows that if this symmetry exists then it must be of the 
form 

V t = fh(V 2 - v 2fi V) + + p^){V x - v lfi V) 

+ Kh{p { ll + p (1) + - U ,_iV) + «/l/l_i(VL 2 - U ,-2^) 

with some indeterminate integration constants «, k. Although the use of this explicit for- 
mula gives no essential advantage in solving the classification problem, it is useful as a final 
check of integrability of the obtained lattices. 

3 Analysis of the integrability conditions 
3.1 First step 

Consider the first pair of integrability conditions (JTj), (jHJ) 

A,(log/) e Im(T- 1), A, (log ft) G Im(T- 1). (13) 
It is easy to obtain the following equations as a corollary: 

/*.-i + fc r (%^ =0, %i + {T^f%i^) =0. (14) 



P f \ P J ' ft 2 ft V ^ 2 

Indeed, the terms containing scalar products Vk,k-z appear only by differentiating with 
respect to a;: 

A,(log/) = h^LD^v^) + ■■■ = ^(A^,-i + ft-i«i,- 2 ) + • • • 

Im( -" 1} ^-v 1 + rf%VU-x + - 



/ J V / 

and the first equation (JHJ) follows. This computation is actually equivalent to applying of 

variational derivative S/5 V30 defined by formula 

*l= » x; t'w, , = 1,2,... 

The use of this notion makes the computations more algorithmic, due to the equality 

5 



C © Im(T - 1) = p| ker 



which is proven along the same lines as in scalar case [2]. 
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Statement 3. The dependence of the coefficients of the lattice on may be one of the 
following: 

Casel. f- aK " l} GKo) 



vi-i + b(v h0 , v ,-i) ' v i _i + 6(«i j0 , uo,-i) ' 

Case 2. f = /(t>i,o, v ,-i), /i = /i(i>i,o, t>o 

Proof. First equation (1T41) implies that f Vl _ 1 / f 2 may depend on i>o,-i only. If /^.j 7^ 
then we come to the Case 1. If f Vl _ ± = then h Vl l = as well, in virtue of the second 
equation (JHJ), and we come to the Case 2. □ 



Conditions (1131) are far from being exhausted by this statement. We will see that in the 
Case 1 they allow to define functions a, b as well. 

3.2 Case 1: j Vx _ x ± 

Notice that in this case the relation ( TTTT) at j = is satisfied with = — log a(t>o,-i). This 
means that conditions ( ITBl) are equivalent to each other and we may consider only the first 
one. Applying of S/S V2 to it is a rather tedious task. The resulting equation is polynomial 
in variables Vk+2,k an d vanishing of the coefficients brings to a certain overdetermined 
system for functions a and b. It is convenient to introduce the auxiliary functions 

1 -v 2 b(u,v)+uv 
y{v) = 2( r , c[u,v) = — 15 

and to denote w = t> 1)0 , v = t>o,-i, «j = U-i,-2- This allows to rewrite the system in a 
relatively compact form as follows: 

c(w, v)(a'(u) - a'(v)) = (a(u)y(u)) u - (a(v)y(v)) v , (16) 

i \i 1 w 1 \i 1 w u(c — y(v)) v(c — y(u)) . . 

a(u)(c + 2/(u))c„ - a(u)(c + y(u))c„ = , c = c{u,v), (17) 

a[u) a{v) 

(c(v , w) + y(v))(2c(u, v) + y(v)) v = (c(u, v) + j/(v))(2c(t;, 10) + y{v)) v . (18) 
At first, we will prove that all solutions of equation ( |T8l) are: 

(i) 2c{u,v) = 2a-y(u)-y{v), 

(ii) c(it, u) = az(u)z(v) + /3, y(v) = ^z 2 {y ) — (3, z' 7^ 

where cn,/3, 7 are arbitrary constants. 

If c(f , w) + y{v) = or c(m, t>) + = then (ITS]) is reduced to the equation 

= (y(u) -y(v))y'(v), 

hence = — (3, c(u,v) = j3, a special case of solution (ii). 

If (c(v, w) + y(v))(c(u, v) + y(v)) 7^ then the variables in (TT51) can be separated: 

(2c(?i,t;)+^)) M = (2c(u,v)+y(v)) v = 

c(u,v)+y(u) c(u,v)+y(v) 
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and as a corollary we obtain c uv = k(u)c v = k(v)c u . The case k = corresponds to 
the solution (i). At k ^ we get c = C(K(u) + K(v)), K' = k, C" = C, whence 
c = az(u)z(v ) + (3, where z' = kz. Moreover, both equations ( fT9l) are reduced to the 
relation 

2z'(v) 

y\v) = ^-±{y{v) + (3) 

and we get (ii) by integration. Now we consider both types of solutions separately and 
come to the following statement. 

Statement 4. The solutions a = a(v), b = b(u,v) of the system ( 1731) -( TlSj) are exhausted, 
up to the scaling a — > const a, by the following list: 

a = v — l/v, b = —uv; (20) 
a 2 — kva + v 2 — 1 = 0, b = a(u)a(v) — uv; (21) 
a = v + £j b = -l; (22) 



a = v + e, b= (u + e)(v + e) fy ^ k ) ~ k ) + k j ~ uv ' ( 23 ) 

a = v + £j b = 1 + e(u + v) + ky/(u + e)(v + e) (24) 

where e = ±1 and k is an arbitrary constant. 

Proof. Solutions of type (i). Applying d u d v to ( TT61) yields 

y'{u)a"{v) = y'(v)a"(u). 

If y' — then scaling allows to set ?/ = 1, a 2 (f) = 1 — v 2 and then (Tl6l) implies that c = 1. 
Equation (fl71) becomes identically true in virtue of these relations and we arrive to solution 
(glD at k = 0. 

If 2/ ^ then a' = fiy + v. The variables in ffTUj) are now separated and we obtain the 
overdetermined ODE system for the functions a = a(v), y = y(v): 

3 

ay' = R{y) = ~-^y 2 + {a^i- v)y + \, a ' = S(y) = fiy + v, ya 2 = l-v 2 . (25) 

Differentiation yields 

a(2yS + R) = -2v, S(2yS + R) + (2S + 2y5 + R)R + 2 = 0. 

The polynomial on y in the l.h.s. of the latter equation must vanish identically since y' 7^ 0. 
This gives the relations \i — 0, \v — — 1 and moreover, the scaling allows to set v — 1. 
Now, the system (|25il is reduced to equations 

ay' = — y — 1, a = v + e, a 2 y = 1 — v 2 . 

It is easy to prove that they are consistent at e 2 = 1, and an intermediate substitution into 
(fTTl) proves that a = 0. The resulting solution is ( 1221 ). 
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Solutions of type (ii). Applying d u d v to ffTBj) yields 

/ ,. . ... z(u)a"(u) z(v)a"(v)\ , . 

a a'fa) - a'fa) + K '-. \ ' - K ' , J = 0. 26 
V z'(v) J v ; 

If a = then c = j3 and variables in equation (fTTj) are separated: 

(/3 - y(u))o(u) (^-J/(u))a(u) 



5. 

This relation turn the equation (fTBl into identity as well. Taking ( fT5l) into account, we 
obtain the equation /5a 2 — 5va + v 2 — 1 = for a(t>). This brings, up to the scaling, to the 
solutions (12U1), (jUD. 

If a 7^ then we set a = 1 without loss of generality. Equation ( |26l) implies a' = n/q+u, 
then the variables in ( fl6l) are separated and we obtain the overdetermined ODE system for 
the functions a = a fa), z = z(v): 

a' = ^ + v, (faz 2 - (3)a)' - — + fiz = X, faz 2 - (3)a 2 = 1 - v 2 . (27) 

Notice that 7 / 0: otherwise —2^(3/z + fiz — [3u = A and since z' 7^ 0, hence fi = 0; but 
then the equations a' = z/, /5a 2 = v 2 — 1 are inconsistent. Therefore, second equation (J2j 
can be rewritten as follows: 

' 1 ( ( , ^ , A + fo/ , 2^/5 

z = - — -71/2 - //(7 + 1) + + — =- 

27a \ z z l 

Now, differentiating of third equation ( 1271) brings, as in the previous case, to a polynomial 
equation for z which must be satisfied identically. This gives equations for the parameters: 

( 7 -l)/5/i = 0, (37-l)(A + /?z/)fz = 0, (7-3)/^ = 0, 4 7 (Az/ + 1) + (7 - 1) V = 0. 

Moreover, substitution into ( fTTj) gives additionally the equations 

(7 + 1)(7 - 3)/5/i = 0, ( 7 2 - 1)(A + = 0, ( 7 2 - l)/i = 0. 

The solutions of the whole system are: 

( /i 2 = l,/5 = 0,A = 0,z/ = 0, 7 = -l), 
fa = 0, v = -1/A, 7 2 = 1), (A* = 0, v = -1/A, /3 = A 2 ). 

The first one is unsuitable since it leads to z' = 0. For the other two we set z/ = 1, A = — 1, 
a = v + e without loss of generality. It is easy to check that (1271) are consistent at e 2 = 1 
and we come to solutions (1231 and (l2"4l . respectively. □ 

It can be proved straightforwardly that conditions ([7J) at j = are fulfilled for each 
solution f[2"U]) - ([2"4"j) . that is there exist quantities a^°\ which turn them into identities. 
It is sufficient to compute only o~(°\ due to the relation a_l = D x (a^) — cr^ where 
&(°> = — log a(f 0i -i). Practically, this computation is based on the "summation by parts" 
algorithm, see e.g. [21 Theorem 1]. After finding one can continue the integrability test 
with the next pair of densities ([9]). It turns out that in all cases except for (1241) the second 
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integrability condition is fulfilled automatically. In the case ( J2IJ) we obtain the restriction 
k 3 — Ak = on the values of parameter. In more details, the density p^ is in this case of 
the form 

(i) f-i i \ I ff-i ( . 

v-x-2 + e v-x-2 + e\ 

~ g (*V^i,o + e + 2£a/^o-i + e) + £ + 2e^/v^~ l 

and it can be proven that 5D x {p 1 ) /5 V20 vanishes if and only if the above constraint holds. 

The computation of and further check of the integrability conditions require the 
considerable efforts. Fortunately, it is possible to avoid these calculations by checking that 
the explicit formula (1121) provide the higher symmetry indeed. This turns out to be true 
for (J20l) -(l23l ) and ( 1241) at k — 0, ±2 (with constants k = —1, R = in all cases) and we 
come, respectively, to the lattices (jVj}~(E0) i n ^ e List [T] below. 



3.3 Case 2: f Vl _ x = 

Computations here are easier, but also more lengthy, since in some subcases we have to 
check up to three integrability conditions (JTj). However, the result of this search is somewhat 
disappointing: it consists of one lattice By this reason we give only schematic account 
of this case. 

Applying S/5 V20 to (fT3|) yields the equations 

^ ( rp f /«D,-n _|_ fvifi \ _|_ fvo,-l _|_ rp-l ( fvi.O 

f\\f)f)f V / 
/I U / & / ^ v ^ 

In turn, differentiating this with respect to t>2,i yields 

(log/k^c-i = °> ( lo g^k, ,vo,-i = f = T(a)b, h = T(c)d 

where a,b,c,d are functions on t>o,-i- Now, the variables in equations ( |28l) are separated 
and we come to relations 

(ab)' c _ {ab)' b _ ( {cd)' b _ {cd)' c _ 
ab a ' afe d cd d , cd a 

with some constants /i, ia If afe + cd 7^ then (a&)' = (cd)' = 0, so that two cases are 
possible, up to the scaling: 

(i) b — p/a, c—ap/p', d——p'/a, p' 7^ 0; 

(ii) a = a/b, d = 1/c. 

In the case (i), applying S/5 V10 to (|T5|) brings to certain overdetermined system for functions 
a, p. It is convenient to analyze this system taking into account some additional information 
(namely, the equation pp" = const {p') 2 ) which can be obtained either from the integrability 
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a{V!-vi fi V ) + a 1 {v 0i - 1 V-V- 1 ) 1 
14 = , a = u„ ! ; (Vi) 

- ^1,0^0-1 U0.-1 

K = ; , a - 2kv _ia + w _ x - 1 = 0; (V 2 ) 

"1,-1 - ^1,0^0-1 + OOi 

^ = K-i + e)(^ + eV) - Ko + e)(K-i + eV) . 

fi-i - 1 

T/ (uo-i + e)(^i + e^)-Ko + e)(Vli + ey) , 

*4 = 7 77 777 7, P=\ ( V4 I 

«i -1 - ^1,0^0 -1 + (^i,o + e)(«o,-i + e)(K+SPi) 




y.= K- 1 + £ )W + ^)-K, + £ )(v-^v) fc = 0±2 . W) 

+ e(v 1>0 + vo-i) + 1 + fcy^ + eV v o,-i + e 

»4 = — F — F > 5 = °> ±L i V 6) 

Vlfl + V ,-l + 

List 1: Integrable lattices, (V, V) = 1, t> m>n = (Kn, Ki), £ = ±1. 



condition ( TTT1) at j = 1 or from the next pair of conservation laws ([2), (EJ). This allows to 
prove that functions a(v),p(v) may be the following: 

1 3 

a = p = -: a — 1, p — v + o; a = v, p = v • . 

v + o 

The check of conservation laws (EJ) , (J9]) for the first solution proves that S must take the 



values ±1,0 and leads to the lattice (\Vf$ , while two other solutions do not pass the test. 

In the case (ii) the first pair of integrability conditions ((7|), (jHJ) is fulfilled for any 
a, b, c. The further analysis proves that conditions ((7j), © are fulfilled if a = 1 and either 
b(v) = c(v) = \/v + 6 or b = c = 1. However, the next conditions ((7|), ffTUl) fail in both 
cases, so that this case turns out to be empty. 

4 The list of integrable lattices 

Theorem 1. If isotropic Volterra type lattice on the sphere (V, V) — 1 satisfies integrability 
conditions lj^)- l[Tu}) then it coincides with one of the lattices from the ListUl up to scaling 
of x. Each lattice from this list possesses at least one higher symmetry of the form {TB) . 



Remark 2. The lattices corresponding to the different signs of e or 5 are equivalent modulo 
flip map V n — > (— l) n V n . The lattice fiV^j at k = ±1 coincides with flVsj ) at k = 0. 



The lattice JT^I ) is the discrete Heisenberg spin chain introduced in [13] . see also [HH 
[T7] where the applications to the discrete geometry were considered and [8J where the 
anisotropic version (see Section [?]) was studied. It can be written (at 5 = 1 and after 
scaling x) as 

Vi+V V + V i 
V = 1 1 29) 
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In this form, the constraint (V, V) — 1 is not necessary for integrability. This lattice and 
its higher symmetry f[T2]) can be written compactly as 



V x = (T - 1){W), V t ={T- l)P w {W 1 - W-i), W = (V + V^Y 1 
by use of the operations 

1 



A' 1 



-A, P A (B) = 2(A,B)A- (A,A)B. 



(A, A) ' 

The variable U satisfies the polynomial lattices 

W x = -Pw{W x - W-!), W t = -P w (P Wl (W 2 + W) - Pw.AW + W- 2 )) 

which are integrable not only in the vector case, but also in more general setting related 
to Jordan triple systems [IS] . 

The lattices flViTHflVsl) are new, up to the author's knowledge. The lattice ( [V3I ) is related 
to flVel ) by composition of difference substitution and reduction. Namely, first we can resolve 
the constraint (V, V) = 1 by use of stereographic projection as explained in Introduction. 
This brings flV^j ) at e — — 1 to the form 



u x = 


\U - 




~U) + \U l -U\ 2 {U -U^) 








and then substitution V = 


u - 


U-i bring 


;s it to the lattice 


v x 


\v 


2 v 1 + \v 1 


2 V \V^\ 2 V + \V\ 2 V^ 




\Vi + v\ 2 


iv + v^i 2 



This is not the same lattice as however it is obvious that both lattices admit the 

reduction on sphere which brings them to the lattice flVgl ). The question on the substitutions 
for the other lattices from the list is so far open. 



5 Associated partial differential equations 

The very general observation due to Levi [TH] is that a higher symmetry of an integrable 
lattice gives rise to some PDE after elimination of the discrete variable n. The lattice itself 
is now interpreted as Backlund transformation for this PDE. The examples of such relation 
can be found in [HI [12] and many other works. In particular, the integrable Volterra lattices 
(jlj) are associated with some systems of nonlinear Schrodinger type. There are known also 
many results on the multifield analogs of NLS-type systems, see e.g. ED [221 [23], 
however their classification is far from being completed. The list of vector Volterra lattices 
provides several new examples of such systems. 

The elimination of the discrete variable is done as follows. The equations (J2J) , ([3]) imply 
the corollaries 

(V x , Vi) = (1 - v 2 )f + (vt _i - v 1>0 v 

(V x , V x ) = (1 - v 2 lfi )f 2 + 2(« 1 ,_ 1 - v lfi v ^)fh + (1 - vl_ x )h 2 . 
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We assume that these equations can be solved with respect to the scalar products f^-i, 
t>o,-i (this is true for all lattices from the List [T]). Then equation can be rewritten in 
the form 

V„ 1 = fV 1 + gV + hV x (30) 
with coefficients depending on the scalar products of vectors Vi,V,V x . Analogously, 

V 2 = fV 1 + gVi + hV. 

Iteration of these formulae allows to express all vectors V n through the vectors U = V\, V 
and their derivatives. As a result, the symmetry (I12p gives rise to a system of the form 

U t = U xx + aU x + f3V x + 7 t7 + SV, 

- (U,U) = (V,V) = 1 (31) 

-V t = V xx + aU x + (3V X + 7*7 + SV, 



with coefficients depending on the scalar products of U, U x , V and V x . The equation ( 1301) 
becomes an explicit Backlund auto-transformation 

U_i = V, V- 1 = fU + gV + hV x 

of this system. Converse is not true: not any integrable system (1311) admits auto-BT of 
such form. Classification problem for this type of equations may be difficult, since even 
the simplest lattices from our list correspond to rather cumbersome systems (13 ip . Few 
instances are given below. In the case flVgj) at 5 = ±1 we come to the system 



_ 2(U x ,V)+46 2V X f (U X ,U X ) 2(U,V X ) \ 

_ 2(U,V X )-45 W _ 2*7 X . ( (V X ,V X ) 2(U X ,V) \ 

* " (u,v) + s Vx (u,v) + 5 + {(u,v) + s + ((u,v) + 5yj [U+0Vh 

while (|V^| ) at 5 = corresponds to the system 

Vl = Utt - mmmi Ux + ( {u „ Ux) + mm v .(™ 



{u,v) 2 V (u,v) J \(u,v) 

. Vt = Vm _ + ((v x , k) - mm v j*> 



{U,V} 2 V (u,v) J \(U,V) 



The lattice (jV^j is associated with the system 



f (U,V x )(U x ,V) _ (u x ,v x -v) \ (U X ,U X ) , 
\{{U,V)+e)* ~ <r/,V)+e ; x (U,V)+e 



+ ju^( 1 + ju^0_\ {£U + V) 

_T/-T/ ■ nf ^l}QS^lYl g^ + gT W rr 

*" M+ V((^^) + £ ) 2 <i/,v) + e ; * + + e * 

+ {V - Vx) (l- {U - V) ) (U + eV) 
(U,V)+e\; (U,V)+e) [U+£V) - 
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6 Presymplectic structure 

The bi-Hamiltonian structure of the scalar Volterra lattice is well known, see e.g. [24J. In 
the vector case the question is more difficult and it requires further investigation. However, 
the following statement shows that all lattices under scrutiny possess at least some uniform 
presymplectic structure. 

Statement 5. Any lattice can ^ e written in presymplectic form 

SV X = ^ + XV, H = = log/K-i.ui.o^o.-i) (32) 

where S is a certain skew- symmetric operator of the form 

S = pT- 1 - pjT - qV-iV y T~ l + q 1 V 1 V T T + r{V x V\ - V-iV?), (33) 

A is Lagrange multiplier corresponding to the constraint (V, V) = 1 and operator UV T acts 
accordingly to the formula UV T (W) = U (V, W) . 



Proof. The equation ( 13211 is equivalent to 

(pT- 1 - Pl T)(fVi +gV + hV-i) - V-iigT- 1 +r)(f + v lfl g + v^h) 
+V 1 (r + qiT^v^f + u o ,-i0 + h) — \V 

_ rp ( /»l,-lT/ I /»0,-1t/ | , /fl,0 T/ . /«0,-l T r . rp— 1 / /«1,0t/ , V 

-T\ CT -V 1 + —V)+—V 1 + —V- 1 + T + 

Equating the coefficients at V, V± 2 yields 

A = pf-i -Pih, p=-fvi,-i/f 2 , pf + Pih = 0. 

The first two equations are just definitions of A and p while the latter one is fulfilled for the 
lattices from the list in virtue of (1141) . Equations for the rest coefficients give the system 
for q and r of the form 

Ar + A x q x = C, Br + B_ x q = D (34) 

where 

A = vi _i/ + Uo.-iP + h, B = f + v lfi g + v^h, 
C = pigi + (log/i/)„ li0 , D = pg^ x - (log . 

Elimination of one of the unknown functions, say r, brings (|34|l to the form 

(T- = BC — AD. 

This means that the system (I34p is solvable if and only if BC — AD G Im(T — 1). Re- 
markably, this condition is equivalent exactly to D x (\ogf) e Im(T — 1), as an easy check 
proves, and therefore it is true for all lattices from the List [U □ 

The concrete expressions for the coefficients q, r may be rather cumbersome (it is clear 
from the proof that they are related somehow with the quantity a^). The answer is very 
simple for the lattice JT3I ): 

11 

p = — , q=- — -, r = 0. (35) 

v -i + e (vo,_i + e) 1 
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The formula (U, SW) = Q(U, W) relates operator 5* with 2-form 



Q = ^{Pn{dV n A dV n -.i) + q n (V n , dV n -x) A (K-i, dV n ) 

n 

+r n (V n+1 ,dV n ) A (V n -i,dV n )) 

where (a A p)(U, W) := (a(U),/3(W)) - (a(W), (3(U)). It is easy to see that this form is 
exact in the case (1551) . namely Q = d^2 n p n {V n , dV n -i). Therefore dfl = 0, that is operator 
S is symplectic indeed. Unfortunately, this is not true in the general case. 

It is also worth to notice that the representation ( J32l) can be replaced with a linear pencil 
by assuming that Hamiltonian is of the form H = p^ + up, where p is some additional 
conserved density depending on v\$ (it does not belong to the sequence (J7|), however it 
turns out that such densities exist for all lattices under consideration). Operator S also 
acquires linear dependence on k, preserving the same structure flHHj) . We bring the explicit 
formulae only for the relatively simple case of lattice §VH : 

(o) i a i 1 1 

P ( ' = log , p = log v i,o, V = -, a= Uo,-i - 



vi-i - Ui,<M),-i a v o,-i 

1 / ,v Ol -1 - ^1,0^0 -l) (^0-2 - ^0-1^-1 -2) 
1 = — + i K - 1) 7 w 

a 2 / ^o,-i A / ^o-i 

<M>0,-1 U -2 

V ' ^1,0 / V U-1,-2 

1 , / ^ ^1,-1 ~ ^1,0^0-1 



aid / «o,-i \ / ^1,0 

Vl,-l 

^i,o / V «b,-i 



Operator S 1 is not symplectic here. We see also that its simplest form corresponds to the 
Hamiltonian p^ + p rather than p(°\ but this may be not so for the other lattices. 



7 Concluding remarks 

The goal of the present paper was to solve some classification problem; such important 
things as difference substitutions, Lax pairs, Backhand transformations, explicit solutions 
and so on have not been considered. These open problems require, probably, more individ- 
ual investigation for each member of the obtained list. From the author's point of view, 
the question on the Hamiltonian properties of the vectorial equations is among the most 
intriguing ones. 

It was mentioned in Introduction that the assumption (iii) can be removed by use of 
stereographic projection. Another interesting setting is related with the variables on the 
cone (V,V) =0 instead of the sphere. At first sight, this constraint may be treated as a 
limiting case, but actually it defines some independent class of equations. In particular, in 
this case the coefficient g is not expressed through /, h and we also have no explicit formula 
like (fl2l) for the symmetry. An interesting example here is the lattice 

V x = —^—(vo-iVi - viflV-i) + b(vx _i, ul, , v -i)V, v njn = 0. 
v i,-i 

It is likely that it satisfies the infinite sequence of integrability conditions (JTj) at arbitrary 
b, but (local) symmetries exist only if b Vl _ x = 0. 
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The other possible generalizations are related with the condition (i). The simplest 
anisotropic lattice is analog of flVg| ) 



/ v : _|_ v y + V i \ 

where K is an arbitrary symmetric operator. This lattice is closely related to many other 
integrable equations, among them Sklyanin lattice and Landau-Lifshitz equation [8]. The 
classification problem in the anisotropic case can be in principle solved along the same lines 
(cf [21 H] in the continuous case), however technically it is much more difficult since coef- 
ficients acquire dependence on the additional variables v m ,n = (V m ,KV n ). It is interesting 
to consider also the asymmetric scalar product (v m , n 7^ f n ,m), however the examples of this 
type are not known at the moment. 

Acknowledgements. This work was supported by the Russian Foundation for Basic 
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